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ABSTRACT 

We present the faintest mid ecliptic latitude survey in the second part of HST archival search for 
outer Solar System bodies. We report the discovery of 28 new trans-Neptunian objects and 1 small 
centaur (R ~ 2km) in the band 5° — 20° off the ecliptic. The inclination distribution of these excited 
objects is consistent with the distribution derived from brighter ecliptic surveys. We suggest that 
the size and inclination distribution should be estimated consistently using suitable surveys with 
calibrated search algorithms and reliable orbital information. 
Subject headings: Kuiper Belt: general - Planets and satellites: formation 



1. INTRODUCTION 

Trans-Neptunian objects (TNOs) represent the left- 
overs of the same planetesimals from which the planets 
in the solar system formed. TNOs offer a unique oppor- 
tunity for testing theories of the growth and collisional 
history of planetes imals and the dynamical evolution of 



the giant pl anets (Kenyon & Bromley 2004 Morbidelli 
et al. 20081. The study of the orbital distribution ot 
i'JN Os has shown the existence of multiple distinct dy- 



namical populations (Levison & Stern 2001 ; Brown 2001 ) 
with differe nt colors (Guibis et al. 2010; Doressoun diram| 
et al.|[2008) and size dist ributions (Bernstein et al.|20TjI 



2008). 
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The population or small (~50 km) TNOs contains mul 
tiple significant clues to understanding the formation of 
the Solar System. There are more faint than bright 
TNOs in the Solar System, which means that faint TNOs 
are a more thorough dynamical tracer, and that some 
subtle clues in the dynamical distribution of TNOs may 
only be revealed by studying this population. Addition- 
ally, Pan & Sari snowed that the size distribution of ob- 
jects at this size, and in particular the size at which the 
size distribution undergoes a change in slope, records the 
collisional history and intrinsic strength of the TNO pop- 
ulation. 

Because of the importance of these small — and there- 
fore faint, with R>25 and in some cases R>27 — TNOs, 
a great deal of effort has been dedicated to searching 
for faint TNOs ([Chiang fc Brown||i999l |Gladman et al. 
2lM||ATIen" et al.|2002||Bernstein et al.|2004[ |Petit et al. 
200m Fraser et al.|2008; Fuentes & Holman|2008| IFraser 



&: Kavelaars|2009[ Fuentes et al.|2009[ |2010p These sur- 
veys have been concentrated near the ecliptic, where the 
sky plane density of objects is largest, since TNOs pass 
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through the ecliptic, regardless of inclination. However 
TNOs with inclination i spend most of their time at eclip- 
tic latitudes ±i. Few deep TNO surveys have been car- 
ried out at ecliptic latitudes greater than a few degrees 
Consequen tly, elaborate debiasing techniques have been 
devel oped ( |Brown|[200T| |Elliot et al.||2005l |Gulbis et aL 
20101 to derive the true inclination distribution ot the 
IN population. There have been no direct measure- 
ments of the inclination distribution of faint TNOs to 
date. 

The Hubble Space Telescope (HST) presents a unique 
opportunity in TNO studies by observing across the en- 
tire sky. Observations made with HST are deep, with a 
single 500 second exposure with the Advanced Camera 
for Surveys (ACS) reaching ~27th magnitude, depending 
on the bandpass. The combination of these two factors 
implies that faint TNOs appear serendipitously in a large 
fraction of all HST images, including fields both on and 
off the ecliptic. The HST archive, therefore, offers the 
opportunity to probe the history of the Solar System by 
measuring the properties of faint TNOs at a wide range 
of ecliptic latitudes. 

We have developed a pipeline that harvests these 
serendipitous TNOs from archival HST/ACS data 
( Fuentes et al.|2010 hereafter F10). In F10, we searched 
ACS data within 5 degrees of the ecliptic and discovered 
14 TNOs, including one binary object. Here we expand 
this search to higher ecliptic latitude, 5-20 degrees. In 
Section 2 we briefly summarize our field selection crite- 
ria and data processing pipeline. Section 3 presents the 
objects discovered in this mid-latitude search and their 
orbital parameter distribution. In Section 4 we use this 
survey to test our current models of inclination distribu- 
tion. 

2. DATA 

We follow the same criteria described in F10 for quick 
identification of objects with TNO-like orbits in data 
taken with ACS/WFC. In this survey we considered 
fields taken at ecliptic latitudes in the range from 5 to 20 
degrees, and for which the total exposure time within a 
pointing was over 1,500 seconds in at least three images. 
The Multimission Archive at STScI (MAST) lists data 
from 1,141 different HST orbits that meet these criteria 
that were available as of June 4th, 2010 (See Fig.[TJ). 
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Fig. 1. — Map of the sky in J2000 coordinates. The region we chose to select our pointings from is shown in red, that in blue is the 
region F10 analyzed. The locations of all data considered are plotted as yellow triangles. The green circles are fields that were not surveyed 
due to problems in the processing or because the sky density yielded too many false positives. Note that there is a concentration of such 
fields near the galactic center (—93.6, —29) where the star density is highest. 



2.1. Pipeline 

In order to allow the implantation of a control popu- 
lation we considered fiat-fielded and undistorted images, 
as well as the filter d ependent PSF models an d field dis- 
tortions explained in |Anderson fc King] ( |2000[ ) . 

The pipeline uses the orbital solution and photometric 
parameters in each image to create and implant the con- 
trol population. After the images have been cleaned of 
cosmic rays (CR) using a Laplacian filter the images are 
ready to be searched for moving objects. 

As the motion on the sky of a TNO varies greatly with 
its orbital elements, we compile a set of elements that 
produce distinct tracks in the images. For each set of 
elements in the set we produce sample PSF tracks that 
are then used to find consistent sources in the images. 
Finally our algorithm links sources among the images in 
the set, finding all those that appear in at least 3 of them. 

Then, a set of images around all plausible combinations 
of detections are created to be scouted by a trained op- 
erator in order to spot obvious contamination by chance 
alignment, elongated objects, poorly subtracted CRs, 
etc. 

2.2. Control Population 

The detection and identification of a TNO in a given 
field of view will depend on the characteristics of the 
observation and the search algorithm. In this project we 
need to characterize our search method to understand 
how efficient it is in finding a moving object in the trans- 
Neptunian region and what variables affect the detection. 

We approach this problem by creating a synthetic pop- 
ulation that is taken from two different orbital distri- 
butions. One is chosen to resemble the semimajor axis 
distribution of the TNO population so that objects will 



look like what we expect are typical TNOs. The second 
distribution tries to sample the whole range of possible 
bound orbits for objects at a distance between 20 and 200 
AU, sampling our ability to detect objects with unusual 
orbits. 

For every original image 200 synthetic objects were in- 
serted. Then the detection process continues with no 
knowledge of which objects are real and which are not. 
Only when real objects are identified and a detection ef- 
ficiency function is constructed is the control population 
position list unveiled. 

3. ANALYSIS 

We run our automated search algorithm consistently 
for all 1141 pointings, of which 98 failed due to too large 
shifts between the images, poor signal to noise, etc. Of 
the remaining 1045 pointings only 986 were visually in- 
spected. A typical field would have less than 500 possible 
detections, so we did not inspect any field that had more 
than 1,000 false positives. The 29 real objects discovered 
and their best fit values for their orbital parameters are 
listed in Table CD 

We use images like those in Figure [2] to evaluate a x 2 
statistic when fitting for the orbital parameters of mov- 
ing objects. This is used in a Markov Chain Monte Carlo 
(MCMC) method to estimate the posterior distribution 
function for the orbital parameters of each found object 
(See Fig. [4] and Fig. [6]) . The only constraints we consid- 
ered on the orbital parameters were a zero velocity along 
the line of sight and a bound Keplerian orbit. 

We usually recovered half of the synthetic objects im- 
planted, yielding a reliable measurement of the bright- 
ness efficiency funciton for every field. The coadded effi- 
ciency function is shown in Fig. [3] This function is well 
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represented by the formula fl e ff = erfc( R 2 ^ 50 ) deg 2 
where Q = 1.83, i? 50 = 26.44 and w = 0.29. 

Among the moving objects discovered there is an ob- 
ject (hst39) whose orbital parameters confidently con- 
strain it to the Centaur region (a < 30AU, Trilling et 
al. in preparation). This object exhibits the least un- 
certainty on its distance, eccentricity and inclination, as 
shown in Table [l] This is due to the improved resolution 
of the motion of HST at the objects' distance. In Fig- 
ure [2] we show the four images taken over the course of 
the pointing coadded. 

4. DISCUSSION 

Of the orbital parameters for which we obtain the best 
constrained quantities are the distance to the object and 
the inclination. The eccentricity is also well constrained 
for nearer objects. 

We present a set of diagrams featuring the inclination 
(Figure E) and eccentricity (Figure [6| as a function of 
distance for the objects in this survey and for the known 
object in the outer solar system (Figure [5] and Figure [7]). 
The information of known objects was taken from the 
JPL Horizons websitcj^] and the distances were computed 
for January 1st, 2011. 

By comparing the shape of the 1 — a confidence re- 
gion in distance, inclination and eccentricity in Figure [4] 
and Figure [6] those quoted in Table [T] we see that nom- 
inal gaussian uncertainties are not representative of the 
actual range of parameters consistent with observations. 
A worse situation occurs in ground based observations, 
with the degrading effect from atmospheric seeing on as- 
trometry. 

It is interesting to note that all inclination solutions 
have a lower bound set by the ecliptic latitude of the 
field where the object was discovered. All but the retro- 
grade solution for hst32 are shown for Fig. Hand Fig. [5] 
While most objects have a relatively well constrained 
distance, data for hst20 is well fit by a wide range of 
distances and inclinations. In Table Q] its nominal dis- 
tance is 79.9 ± 15.3 AU which yields H = 6.9, the lowest 
absolute magnitude among our objects. Comparing the 
probablity distribution for hst20's orbital parameters in 
Figure |4| and Fi gure \6\ we see that the lowest inclina- 
tion corresponds to the highest eccentricity. This solu- 
tion puts hst20 much closer and in a region with already 
known objects as can be seen in Figure [5] and Figure [7j 

Some trends can be readily observed, like the high con- 
centration of objects at around ~ 42 AU in Figure [4] and 
Figure [6| The eccentricity is not well contrained except 
for the Tew objects closer to the observer. Nevertheless, 
the data constrains tight relationships between the or- 
bital parameters. 

In Figure[7]we see a stream of objects whose eccentric- 
ity decreases with distance. These correspond to Pluti- 
nos in 3:2 resonance with Neptune. There is compelling 
evidence for the detection of some Plutinos in our sam- 
ple by looking at Figure [6] where hstl5, hstl8, hst25, 
hst27, and hst31 have distances and eccentricities con- 
sistent with being in the 3:2 resonance. 

In the inclination vs. distance diagrams we excluded 
the one retrograde solution allowed by our data (See 
hst32 in Table]!]). That solution is in an area of the pa- 

2 http://ssd.jpl.nasa.gov/ 



rameter space completely void of known objects, while 
the closer prograde solution (See Table [lj is well accom- 
panied by dynamically hot TNOs. 

4.1. Inclination Distribution 

Our deep, mid-latitude survey allows us to measure the 
inclination distribution of small TNOs. We assume there 
are only two classes of objects in our sample of TNOs, 
dynamically cold and hot (those with low and high incli- 
nations respectively) . We consider a gaussian inclination 
distribution f(i), and a latitud i nal di stribu tion L(j3) like 
the on es described in Brown (2001) and Gulbis et al. 

polo I. 



f j (i)=sm(i) exp(-^|) , 
V sin 1 — sin p 

h 



(1) 

(2) 
(3) 



where j represents the hot or cold population and Kj 
is a normalization constant such that Lj(0) — 1. The 
integrand pj{i\(3) is the expected inclination distribution 
from a single observation made at a latitude /3. 

We consider a population of objects with distances like 
those of classical TNOs and inclinations and eccentric- 
ities drawn from uniform distributions to compute the 
total effective survey area. For two or more pointings of 
the same region of the sky that were observed within a 
couple of days we calculate the total effective area and 
detection efficiency depending on the number of objects 
in our test population that appear in multiple pointings. 
The result is a set of independent effective observations 
with different detection efficiencies and observed areas. 
A total of 777 such observations were identified and vi- 
sually inspected from the 986 pointings searched. 

We then compute the expected number of objects for 
our survey as a function of inclination considering both 
the luminosity function and inclination distribution of 
hot and cold TNOs. We compare this to that of the 
TNOs discovered in our survey: 



N „ 
fe=0 3 J -' 



o-j(R)dR(4) 



where o-j(R) is the number density of hot or cold objects 
at the ecliptic (/3 = 0), with the values presented in F10, 
and rjk is the detection efficiency function for pointing k. 

The observed inclination distribution in Figure [8] is the 
addition of the inclination probability distribution for all 
28 TNOs in Table [Tj We fit for the width of the hot pop- 
ulation inclination distribution ah by comparing the ex- 
pected and observed probability distributions using the 
Kolmogorov-Smirnov statistic. The best fit parameter is 
given by ah — 16.5l 3 ; 5 degrees. The effect of consider- 
ing the cold population is marginal, with less than one 
cold object expe cted in our survey w he ther a c = 2° or 3° 



(as computed in Gulbis et al. (2010) or Brown (2001 ) for 
larger TNOs). This result indicates the TJNO inclination 
distribution is independent of size. 
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Fig. 2. — Coadded image of the Centaur hst39, d= 12.9 ± 0.3 AU, over a single HST visit showing the apparent motion of the object. 
The angular size of the Earth at that distance is shown as a red circle. The white bar indicates the relative size of the off ecliptic motion 
of HST during the visit. The four individual frames are shown in sequence. 

There is a caveat in the the expected distribution of 
objects presented in Figure [8] The expected number of 
hot objects was 50% larger than those in our survey. 
To better compare expected with observed inclination 
distributions we adjusted the density per unit area for 
hot objects (£23) accordingly, from 0.9 to 0.5. This only 
affects the number of expected hot objects, and not the 
shape of the distribution. 



inclination distribution = 16.5^3'j? degrees. Changing 
this parameter to the 1-a confidence limits gives a signifi- 
cantly poorer representation of the observed distri bution. 
The fa ct this result is so close to that found by |Brown 
(2001 1 is an indication that both large and small objects 
have the same inclination distribution. Comparison to 
th e Deep Eclipt i c Surv ey (DES) inclination distributions 



Gulbis et al. (2010) are difficult since we do not have 



5. CONCLUSIONS 

We have conducted a deep search sensitive to TNOs 
as faint as R~27 between ecliptic latitudes 5° — 20°. 
These are overwhelmingly dynamically excited objects 
(only one cold object is expected to appear in our sam- 
ple), in contrast to surveys performed in the ecliptic 
where most of the found objects are not dynamically 
excited. This survey presents a unique opportunity to 
compare the characteristics of dynamically excited ob- 
jects to those extrapolated from surveys carried out near 
the ecliptic. 

The inclination distribution of TNOs has been modeled 
primarily wit h objects discovered near the ecliptic plane 
Brown| pOOl] ) ; |Gulbis etaLl (12010 b ■ The double gaussian 
model introduced by Brown 1 20U1 ) recognized the exis- 
tence of two different populations in the classical TNOs, 
one dynamically excited and the other with seemingly 
unperturbed near- circular orbits. Further differences in 



surface properties jDorcssoundir am et al.j (2008) a nd lu- 
minos it y function [B ernstein et al. 1 2004[) ; |Fuentes et al. 



(2010); Fraser et al.| ( |2010p have been reported in the 
literature, hinting on an underlying difference in size dis- 
tribution and origin. 
From our survey we can constrain the width of the hot 



enough orbital information to identify our objects in each 
of the DES categories. 

There are about 500 TNOs discovered in deep, well 
callibratcd surveys. This is enough to fit the luminosity 
function and inclination distribution in a self consistent 
manner for the hot, cold, and Plutino populations. The 
main difficulty for this is the poor constraint on the in- 
clination that is obtained from deep ground based sur- 
veys. This can lead to confusion and could explain the 
discrepancy in the overall number density between those 
discovered from the ground (that carry most of the sta- 
tistical weight in the F10 luminosity function) and our 
survey. 

The single centaur found at d = 13AU in this survey 
(hst39) has an estimated radius of 2km if an albedo typ - 
ical of TNOs p~7% is assumed ( |Stansberry et al.|2008[ ). 
This is one of the smallest object found in the outskirts 
of the Solar System, close i n size to the 500-m sized TNO 
discovered by occultations ( Schlichting et al."|20 09 ) . This 
discovery places a measurement of the Centaur size dis- 
tribution at a size that is comparable to the known mem- 
bers of the Jupiter family of comets ( JFCs, Trilling et al. 
in preparation). 
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Fig. 3. — Total efficiency function for all pointings considered in this survey. It is well represented by the formula fl e ff = 
n erfc( fl ~^ 50 ) deg 2 where Q = 1.83, R 50 = 26.43 and w = 0.29. 

Support for program 11778 was provided by NASA for Research in Astronomy, Inc., under NASA contract 

through a grant from the Space Telescope Science Insti- NAS 5-26555. 
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Fig. 5. — Distance on Jan 1st. 2011 versus inclination. Outer solar system objects (a>30AU, black squares), centaurs (green crosses) 
and Jupiter trojans (yellow dots) as denned in the JPL lists are shown for comparison. The retrograde solution for hst32 is not shown. 
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Fig. 6. — Distance versus eccentricity la probability contour for each one of the objects in Table ^ The confidence limits reported in 
that table are shown in red. 




100 



4AU] 

Fig. 7. — Distance on Jan 1st, 2011 versus inclination. Outer solar system objects (a>30AU, black squares), centaurs (green crosses) 
and Jupiter trojans (yellow dots) as defined in the JPL lists are shown for comparison. 
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Fig. 8. — The inclination probability distribution for all objects found in this survey is shown as a continuous black line. The expected 
probability distribution for different models consider the luminosity function in F10 and <r c = 3°. The best fit for = 16.5° is shown as a 
blue solid line. The limits for the 1 — a confidence level are also shown as a dashed green line and a dashed-dotted red line (tr^ = 13°, 21° 
respectively). 
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Note. — All objects found in this work are shown with their photometric and astromctric properties. Positions given for the fir st detections. 
The barycentric distance d and inclination i were estimated from a MCMC with a parameterization given by the Orbfit code (Bernstein Sz 
Khushalani 2000). Though some objects were discovered in the same field, the epoch of the observations is different. The Solar System magnitude 
H — V + 51ogd A, a function of the V magnitude d and the distance to the observer A, is computed assuming the phase angle is small and that 
the V-R color for all objects is 0.6. The conversion between HST filters and the Johnson system arc detailed in F10. 
a When prograde and retrograde solutions are possible we report both peaks. 



